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A measurement of the ZZ production cross section in proton-proton collisions at y/s = 7 TeV using 
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corresponding to an integrated luminosity of 4.6 fb _1 collected in 201 1, events are selected that are 
consistent either with two Z bosons decaying to electrons or muons or with one Z boson decaying 
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Standard Model prediction of 5.891^8 pb calculated at next-to-leading order in QCD. The normalized 
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event yield as a function of the transverse momentum of the leading Z boson is used to set limits on 
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Measurement of ZZ production in pp collisions at 

= 7 TeV and limits on anomalous ZZZ and ZZj 
couplings with the ATLAS detector 



Abstract: A measurement of the ZZ production cross section in proton-proton collisions at y/s = 
7 TeV using data recorded by the ATLAS experiment at the Large Hadron Collider is presented. In a 
data sample corresponding to an integrated luminosity of 4.6 fb _1 collected in 2011, events are selected 
that are consistent either with two Z bosons decaying to electrons or muons or with one Z boson 
decaying to electrons or muons and a second Z boson decaying to neutrinos. The ZZ^ — > 
and ZZ — > i^lrvv cross sections are measured in restricted phase-space regions. These results are 
then used to derive the total cross section for ZZ events produced with both Z bosons in the mass 
range 66 to 116 GcV, = 6.7 ± 0.7 (stat.) t A 3 (syst.) ± 0.3 (lumi.) pb, which is consistent 
with the Standard Model prediction of pb calculated at next-todeading order in QCD. The 

normalized differential cross sections in bins of various kinematic variables are presented. Finally, the 
differential event yield as a function of the transverse momentum of the leading Z boson is used to set 
limits on anomalous neutral triple gauge boson couplings in ZZ production. 
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1 Introduction 

The production of pairs of Z bosons at the Large Hadron Collider (LHC) provides an excellent op- 
portunity to test the predictions of the electroweak sector of the Standard Model (SM) at the TeV 
energy scale. In the SM, Z boson pairs can be produced via non-resonant processes or in the decay 
of Higgs bosons. Deviations from SM expectations for the total or differential ZZ production cross 
sections could be indicative of the production of new resonances decaying to Z bosons or other non-SM 
contributions. 
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Non-resonant ZZ production proceeds at leading order (LO) via t- and u-channel quark-antiquark 
interactions, while about 6% of the production proceeds via gluon fusion. The ZZZ and ZZj neutral 
triple gauge boson couplings (nTGCs) are absent in the SM, hence there is no contribution from s- 
channel qq annihilation at tree level. These different production processes are shown in figure 1. At 
the one-loop level, nTGCs generated by fermion triangles have a magnitude of the order of 10~ 4 [1]. 
Many models of physics beyond the Standard Model predict values of nTGCs at the level of 10 -4 to 
10~ 3 [2]. The primary signatures of non-zero nTGCs are an increase in the ZZ cross section at high ZZ 
invariant mass and high transverse momentum of the Z bosons [3] . ZZ production has been studied 
in e + e~ collisions at LEP [4-8], in pp collisions at the Tevatron [9-12] and recently in pp collisions at 
the LHC [13, 14]. No deviation of the measured total cross section from the SM expectation has been 
observed, and limits on anomalous nTGCs have been set [8, 9, 13, 14]. In searching for the SM Higgs 
boson, the ATLAS and CMS collaborations observed recently a neutral boson resonance with a mass 
around 126 GeV [15-17]. A SM Higgs boson with that mass can decay to two Z bosons only when 
at least one Z boson is off-shell, and even in this case, the contribution is less than 3%. Searches for 
high-mass non-SM ZZ resonances have not resulted in any excess above the SM expectations [18]. 




(d) (c) 

Figure 1. Leading order Feynman diagrams for ZZ production through the qq and gg initial state at hadron 
colliders. The s-channel diagram, (c), contains the ZZZ and ZZy neutral TGC vertices which do not exist in 
the SM. 

This paper presents a measurement of ZZ production 1 in proton-proton collisions at a centre-of- 
mass energy y/s = 7 TeV using 4.6 fb _1 of integrated luminosity collected by the ATLAS detector at 
the LHC. ZZ events are selected in two channels 2 : £ + £~£' + £'~ and l + l~vv. Two selections are used 

1 Throughout this paper Z should be taken to mean Z/f* when referring to decays to charged leptons, and just Z 
when referring to decays to neutrinos. 

2 £ represents either electrons or muons. I and £' are used to denote leptons from a different Z parent, but not 
necessarily of different flavour. Decay modes mentioned with the use of £ indicate the sum of the decay modes with 
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in the four-charged-lepton channel: an on-shell ZZ selection denoted by ZZ — > £ + £~£' + £'~ where 
both Z bosons are required to be within the mass range 66-116 GeV 3 and a selection which includes 
an off-shell Z boson denoted by ZZ* —> £ + £~£' + £'~ where one Z boson is required to be within this 
mass range and the other can be off-shell and have any mass above 20 GeV. In the £ + £~vv channel, 
the vv system is expected to be produced by an off-shell Z boson in 2.6% of the events. Since this 
fraction is small and only one event selection is used for this channel, it is referred to as ZZ — > £ + £~ vv 
throughout the paper. The ZZ^ — > £ + £~ £' + £'~ channel has an excellent signal-to-background ratio, 
but it has a branching fraction six times lower than the ZZ — > £ + £~vv channel; the latter has higher 
background contributions with an expected signal-to-background ratio around one (after applying the 
event selections described below). This paper presents the total ZZ production cross section, the 
fiducial cross section in a restricted phase space for each decay channel (integrated, and as a function 
of kinematic parameters for the ZZ selections) and limits on anomalous nTGCs using the observed 
ZZ event yields as a function of the transverse momentum of the leading Z boson 4 . The results 
presented in this paper supersede the previously published results [13] which were derived with the 
first 1.02 fb _1 of the dataset used here, only with the ZZ — > £ + £~ £' + £'~ decay channel and with the 
use of the total ZZ event count for the derivation of the limits on anomalous nTGCs. 

The total cross section for non-resonant ZZ production is predicted at next-to-leading order 
(NLO) in QCD to be 6.18^;^ pb, where the quoted theoretical uncertainties result from varying 
the factorization and renormalization scales simultaneously by a factor of two whilst using the full 
CT10 parton distribution function (PDF) error set [19]. The cross section is calculated in the on-shell 
(zero-width) approximation using MCFM [20] with CT10; it includes a 5.8% contribution from gluon 
fusion. When the natural width of the Z boson is used and both Z bosons are required to be within 
the Z mass window, the NLO cross section is predicted to be 5. 89j^i8 pb. The cross sections given 
here are calculated at a renormalization and factorization scale equal to half the mass of the diboson 
system. The total cross section using the zero-width approximation was previously measured to be 
8.5±|;| (stat.) ta.t (syst.) ± 0.3 (lumi.) pb [13]. 

This paper is organized as follows: an overview of the ATLAS detector, data, signal and back- 
ground Monte Carlo (MC) samples used for this analysis is given in section 2; section 3 describes 
the selection of the physics objects; section 4 describes the fiducial phase space of the measurement, 
the corresponding ZZ cross section definition and the acceptances of the event selection and fiducial 
phase space; section 5 explains how the backgrounds to the £ + £~£' + £'~ and £ + £~vv final states are 
estimated with a combination of simulation and data-driven techniques; section 6 presents the results: 
cross section, differential cross sections and nTGC limits; finally, a summary of the main results is 
given in section 7. 

2 The ATLAS detector and data sample 

The ATLAS detector [21] is a multipurpose particle detector with a cylindrical geometry. It consists 
of inner tracking devices surrounded by a superconducting solenoid, electromagnetic and hadronic 
calorimeters and a muon spectrometer with a toroidal magnetic field. The inner detector, in com- 
bination with the 2 T field from the solenoid, provides precision tracking of charged particles in the 

specific lepton flavours. 

3 Throughout this paper, the 66—116 GeV mass range is referred to as the Z mass window. 

4 Leading Z refers to the Z with the higher transverse momentum in ZZ — ¥ 1+ 1~ l ,J <~ l'~ decays or to the Z boson 
decaying to a charged lepton pair in ZZ — > £ + £~vu decays. 
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pseudorapidity range \r)\ < 2.5°. It consists of a silicon pixel detector, a silicon microstrip detector and 
a straw tube tracker that also provides transition radiation measurements for electron identification in 
the pseudorapidity range \rj\ < 2.0. The calorimeter system covers the pseudorapidity range \r}\ < 4.9. 
The electromagnetic calorimeter uses liquid argon (LAr) as the active material with lead as an ab- 
sorber (|?7| < 3.2). It identifies electromagnetic showers and measures their energy and position; in 
the region < 2.5 it is finely segmented and provides electron identification in conjunction with the 
inner detector which covers the same 77 region. Hadronic showers are measured in the central rapidity 
range (|?7| < 1.7) by scintillator tiles with iron absorber, while in the end-cap region (1.5 < \r)\ < 3.2) 
a LAr calorimeter with a copper absorber is used. In the forward region (3.2 < \rj\ < 4.9) a LAr 
calorimeter with a copper absorber for the first layer and tungsten for the last two layers is used 
for both electromagnetic and hadronic showers. All calorimeters are used to measure jets. The 
muon spectrometer surrounds the calorimeters; it consists of superconducting air-core toroid magnets, 
high-precision tracking chambers which provide muon identification and tracking measurement in the 
pseudorapidity range \rj\ < 2.7, and separate trigger chambers covering |7y| < 2.4. 

A three-level trigger system selects events to be recorded for offline analysis. The events used 
in this analysis were selected with single-lepton triggers with nominal transverse momentum (pt) 
thresholds of 20 or 22 GeV (depending on the instantaneous luminosity of the LHC) for electrons and 
18 GeV for muons. The efficiencies of the single-lepton triggers have been determined as a function 
of lepton pseudorapidity and transverse momentum using large samples of Z — > £ + £~ events. The 
trigger efficiencies for events passing the offline selection described below are all greater than 98%. 

The measurements presented here uses the full data sample of proton-proton collisions at i/s = 
7 TeV recorded in 2011. After data quality requirements, the total integrated luminosity used in the 
analysis is 4.6 fb _1 with an uncertainty of 3.9% [22]. 

2.1 Simulated data samples 

Monte Carlo simulated samples cross-checked with data are used to calculate several quantities used in 
this measurement, including acceptance, efficiency and some of the background to the ZZ — > £ + £~i>D 
decay channel. The NLO generator PowhegBox [23, 24] with the CT10 PDF set, interfaced to 
Pythia [25], is used to model the signal for both channels. The LO multi-leg generator Sherpa [26] 
with the CTEQ6L1 PDF set [27] in comparison with PowhegBox is used to evaluate systematic 
uncertainties. The contribution from gg —> ZZ is modelled by the GG2zz generator [28] interfaced 
to Herwig [29] to model parton showers and to Jimmy [30] for multiparton interactions. In each 
case, the simulation includes the interference terms between the Z and 7* diagrams. For both the 
£ + £~£' + £'~ and £ + £~vv final states, MCFM is used to calculate theoretical uncertainties, and Sherpa 
is used for the generation of signal samples with neutral triple gauge couplings. 

The LO generator Alpgen [31] with CTEQ6L1 PDFs is used to simulate Z+jets, VF+jets, and 
FT 7 background events with Jimmy used for multiparton interactions and Herwig for parton showers. 
The NLO generator MC@NLO [32] with CT10 PDFs is used to model tt background processes as well 
as WW production. The single-top Wt process is modelled with AcerMC [33] with the MSTW2008 
PDFs [34]. The LO generator Herwig with MSTW2008 PDFs is used to model WZ production. The 
LO generator Madgraph [35] with CTEQ6L1 PDFs is also used to model Zj and Wj* events, where 
Pythia is used for hadronization and showering. 

5 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the centre of the 
detector and the z-axis along the beam direction. The x-axis points from the interaction point to the centre of the LHC 
ring, and the y-axis points upwards. Cylindrical coordinates (r,<j>) are used in the transverse plane, <j> being the azimuthal 
angle around the beam direction. The pseudorapidity r\ is defined in terms of the polar angle 9 as r) = — lntan(6*/2). 
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The detector response is simulated [36] with a program based on Geant4 [37]. Additional inelastic 
pp events are included in the simulation, distributed so as to reproduce the number of collisions per 
bunch-crossing in the data. The detector response to interactions in the out-of-time bunches from pile- 
up is also modelled in the simulation. The results of the simulation are corrected with scale factors 
determined by comparing efficiencies observed in data to those in the simulated events, and the lepton 
momentum scale and resolution are finely adjusted to match the observed dilepton spectra in Z — > ££ 
events using a sample of Z bosons. 

3 Event reconstruction and selection 

Events are required to contain a primary vertex formed from at least three associated tracks with 
p T > 400 MeV. 

3.1 Leptons, jets and missing energy 
3.1.1 Common lepton selection 

Muons are identified by matching tracks (or track segments) reconstructed in the muon spectrometer to 
tracks reconstructed in the inner detector [38] . The momenta of these combined muons are calculated 
by combining the information from the two systems and correcting for the energy deposited in the 
calorimeters. The analyses of both decay channels use muons which have full tracks reconstructed in 
the muon spectrometer with p T > 20 GeV and \rj\ < 2.5. The ZZ& -> £+£-£'+£'- channel recovers 
additional ZZ acceptance with minimal additional background using a lower threshold of pt > 7 GeV 
and by accepting muons with segments reconstructed in the muon spectrometer (in this latter case, the 
muon spectrometer is used to identify the track as a muon, but its momentum is measured using the 
inner detector; for the purposes of the discussion below, these muons are also referred to as combined 
muons) . 

Electrons are reconstructed from an energy cluster in the electromagnetic calorimeter matched to a 
track in the inner detector [38] ; the transverse momentum is computed from the calorimeter energy and 
the direction from the track parameters measured in the inner detector. The electron track parameters 
are corrected for bremsstrahlung energy loss using the Gaussian-sum filter algorithm [39]. Electron 
candidates in the ZZ^ — > £ + £~£' + £'~ (ZZ — > £ + £~vv) channel are required to have longitudinal and 
transverse shower profiles consistent with those expected from electromagnetic showers, by satisfying 
the loose (medium) identification criteria described in ref. [40] reoptimized for the 2011 data-taking 
conditions. They are also required to have a transverse momentum of at least 7 (20) GeV and a 
pseudorapidity of \rj\ < 2.47. 

In order to reject non-prompt leptons from the decay of heavy quarks and fake electrons from 
misidentified jets (charged hadrons or photon conversions), all selected leptons must satisfy isolation 
requirements based on calorimetric and tracking information and must be consistent with originating 
from the primary vertex. For the calorimetric isolation the scalar sum of the transverse energies, YjE^, 
of calorimeter deposits inside a cone around the lepton, corrected to remove the energy from the lepton 
and from additional interactions (pile-up), is formed. In the ZZ^ -t £ + £~£' + £'~ (ZZ -t £ + l~vv) 

channel, the EE't inside a cone of size AR = J (A(f>) 2 + (Ar/) 2 = 0.2 (0.3) around the lepton is 
required to be no more than 30% (15%) of the lepton pr- For the track isolation, the scalar sum of 
the transverse momenta, £pr> of inner detector tracks inside a cone of size AR — 0.2 (0.3) around the 
lepton is required to be no more than 15% of the lepton p^. The wider cone size, in conjunction with 
the same or tighter requirements on the fraction of extra activity allowed in the cone, corresponds 
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to more stringent isolation requirements applied to the ZZ — ¥ £ + £~vv channel compared to the 
ZZ^*> — > £ + £~ £' + £'~ channel. This reflects the need to reduce the much higher reducible background 
(predominantly from Z+jets, tt and WW). To ensure that the lepton originates from the primary 
vertex, its longitudinal impact parameter \zq\ is required to be less than 2 mm, and its transverse 
impact parameter significance (the transverse impact parameter divided by its error), Ido/a^J, is 
required to be less than 3.5 (6) for muons (electrons). Electrons have a worse impact parameter 
resolution than muons due to bremsstrahlung. 

Since muons can radiate photons which may then convert to electron-positron pairs, electron 
candidates within AR = 0.1 of any selected muon are not considered. If two electron candidates are 
within AR = 0.1 of each other, the one with the lower pr is removed. 

3.1.2 Extended-lepton selection 

Two additional categories of muons are considered for the ZZ^*' — > £ + £~ £' + £'~ channel: forward 
spectrometer muons with 2.5 < \rj\ < 2.7 (in a region outside the nominal coverage of the inner 
detector) and calorimeter-tagged muons with \rj\ < 0.1 (where there is a limited geometric coverage 
in the muon spectrometer). Forward spectrometer muons are required to have a full track that is 
reconstructed in the muon spectrometer; if these muons are also measured in the inner detector, 
their momentum is measured using the combined information; otherwise, only the muon spectrometer 
information is used. In either case, such muons are required to have pt > 10 GeV and the Y>Et of 
calorimeter deposits inside a cone of size AR = 0.2 around the muon is required to be no more than 
15% of the muon pi>, while no requirement is made on Ep-r- The same impact parameter requirements 
as for the combined muons are imposed for the forward muons measured in the inner detector; no 
such requirement is imposed on those measured in the muon spectrometer only. Calorimeter-tagged 
muons are reconstructed from calorimeter energy deposits consistent with a muon which are matched 
to an inner detector track with > 20 GeV and are required to satisfy the same impact parameter 
and isolation criteria as for the combined muons. 

The ZZ^*> — > £ + £~£' + £'~ channel also uses calorimeter-only electrons with 2.5 < \rj\ < 3.16 and 
Pt > 20 GeV passing the tight identification requirements [40] for this forward r\ region, where only the 
longitudinal and transverse shower profiles in the calorimeters are used for their identification. Their 
transverse momentum is computed from the calorimeter energy and the electron direction, where the 
electron direction is computed using the primary vertex position and the shower barycentre position 
in the calorimeter. Being identified outside the acceptance of the inner detector, no impact parameter 
requirements can be applied to these calorimeter-only electron candidates, and their charge is not 
measured. Since only one such electron is allowed in the event, and since all other leptons have their 
charge measured, the calorimeter-only electron is assigned the charge needed to have two pairs of same- 
flavour opposite-sign leptons in the event. The requirements described above constrain the additional 
background introduced by the inclusion of calorimctcr-only electrons, and no isolation requirements 
are imposed on such electrons. 

The use of the extended-lepton selection increases the ZZ -> £ + £r£' + £'~ and ZZ* -> e + £~£'+t~ 
acceptance by about 6% from the forward spectrometer muons, 4% from the calorimeter-tagged muons 
and 6% from the forward electrons. The expected background is kept small by requiring each event to 
have at most one lepton from each extended-lepton category, and each such lepton to be paired with 
a non-extended lepton. 



- 6 - 



3.1.3 Jets and missing transverse momentum 

For the ZZ — > £ + £~vv selection, events which contain at least one well-reconstructed jet are vetoed 
to reduce background from top-quark production. Jets are reconstructed from topological clusters of 
energy in the calorimeter [41] using the anti-/c t algorithm [42] with radius parameter R = 0.4. The 
measured jet energy is corrected for detector inhomogeneities and for the non-compensating nature of 
the calorimeter using px~ and 77-dependent correction factors based on Monte Carlo simulations with 
adjustments from in-situ measurements [43, 44]. Jets are required to have pt > 25 GeV and < 4.5. 
In order to minimize the impact of jets from pile-up at high luminosity, the jet vertex fraction is 
required to be at least 0.75; the jet vertex fraction is defined as the sum of the pt of tracks associated 
to the jet and originating from the primary vertex, divided by the sum of the px of all the tracks 
associated to the jet. If a reconstructed jet and a lepton lie within AR = 0.3 of each other, the jet is 
not considered in the analysis. 

The missing transverse momentum _E™ ISS is the imbalance of transverse momentum in the event. 
A large imbalance in the transverse momentum is a signature of the ZZ — > £ + £~vD decay channel. The 
two-dimensional E™ lss vector is determined from the negative vectorial sum of reconstructed electron, 
muon and jet momenta together with calorimeter cells not associated to any object [45]. Calorimeter 
cells are calibrated to the jet energy scale if they are associated with a jet and to the electromagnetic 
energy scale otherwise. Using calorimeter timing and shower shape information, events that contain 
jets with px > 20 GeV and not originating from proton-proton collisions but from e.g. calorimeter 
signals due to noisy cells are rejected. 

3.2 ZZW £+£-£'+£'- selection 

ZZ^*' — > £ + £~£' + £'~ events are characterized by four high-px, isolated electrons or muons, in three 
channels: e + e~e + e~, /i + p~/i + /i~ and e + e~ [i + p~ ■ Selected events are required to have exactly four 
leptons and to have passed at least a single-muon or single-electron trigger. Each combination of 
lepton pairs is required to satisfy AR(£\,£-2) > 0.2, where £\ and £2 are used hereafter to denote a 
pair of distinct leptons, independent of their Z parent assignment, flavour and charge. To ensure 
high and well-measured trigger efficiency, at least one lepton must have px > 20 GeV (25 GeV) for 
the offline muon (electron) and be matched to a muon (electron) reconstructed online by the trigger 
system within AR = 0.1 (0.15). 

Same-flavour, oppositely-charged lepton pairs are combined to form Z candidates. An event must 
contain two such pairs. In the e + e~e + e - and fi + p~p + (i~ channels, ambiguities are resolved by 
choosing the combination which results in the smaller value of the sum of \mi+£- — mz\ for the two 
pairs, where nip+e- is the mass of the dilepton system and mz is the mass of the Z boson [46]. Figure 2 
shows the correlation between the invariant mass of the leading (higher px) and the sub-leading (lower 
Pt) lepton pair. The events cluster in the region where both masses are around mz- At least one 
lepton pair is required to have invariant mass within the Z mass window, 66 < m^+g- < 116 GeV. If 
the second lepton pair satisfies this as well, the event is classified as a ZZ event; if the second pair 
satisfies m e + e - > 20 GeV, the event is classified as a ZZ* event. 

With the selection described here, 84 ZZ* — > £ + £~£' + £'~ candidates are observed, out of which 
66 are classified as ZZ -> £+£-£'+£'- candidates. From the 84 (66) ZZ* £+£-£'+£'- (ZZ 
£ + £~£' + £'~) candidates, 8 (7) candidates contain extended leptons. 
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Figure 2. The mass of the leading lepton pair versus the mass of the sub-leading lepton pair. The events 
observed in the data are shown as solid circles and the ZZ^ — > £ + £~£' + £'~ signal prediction from simulation 
as boxes. The size of each box is proportional to the number of events in each bin. The region enclosed by 
the solid (dashed) lines indicates the signal region defined by the requirements on the lepton-pair masses for 
ZZ (ZZ*) events, as defined in the text. 



3.3 ZZ 1+l-vD selection 

ZZ — > £ + £~vD events are characterized by large missing transverse momentum and two high-px, 
isolated electrons or muons. Selected events are required to have exactly two leptons of the same flavour 
with 76 < mi+£- < 106 GeV and to have passed at least a single-muon or a single-electron trigger. 
The mass window is chosen to be tighter than the mass window used for the ZZ^> -> £+£-£'+£'- 
channel in order to reduce the background from tt and WW. The lepton pair is required to have 
AR(£ + , l~) > 0.3. This requirement reflects the choice of the isolation cone for the leptons. The same 
trigger matching requirement as in the ZZ^ -> £+£-£'+£'- channel is used. 

The ZZ — > £ + £~vv decay channel analysis makes use of several selections to reduce background. 
The largest background after the mass window requirement consists of ^+jets events, which are 
associated with non-zero missing transverse momentum when the i?™ lss is mismeasured or when a 
6-quark decays to leptons and neutrinos inside of a jet. Since the Z bosons tend to be produced 
back-to-back, the axial- E™ ss (defined as the projection of the E™ lss along the direction opposite to 
the Z — > t + £~ candidate in the transverse plane) is a powerful variable to distinguish ZZ — > £ + £~vv 
decays from Z+jets. The axial- E™ lss is given by —E™ lss ■ p z /p^, where p z is the magnitude of the 
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transverse momentum of the Z candidate. Similarly, the fractional px difference, |_E™ 1SS — pfl/P? 
is a good variable to distinguish the two. The axi&\-E™ lss and fractional px difference are shown in 
figure 3. In order to reduce Z+jets background, the &xial-E™ lss must be greater than 75 GeV, and 
the fractional px difference must be less than 0.4. To reduce background from top-quark production, 
events which contain at least one reconstructed jet with px > 25 GeV and \rj\ < 4.5 are rejected. 

To reduce background from WZ production, events with a third lepton (electron or muon) with 
px greater than 10 GeV are rejected. The shape of the jet multiplicity distribution is well modelled in 
Monte Carlo simulation as shown in figure 4 for the ZZ —> £ + £r £' + £'~~ and ZZ —> £ + £~vv selections, 
however, there is an overall excess of about 20% in the ZZ — > £ + £r£' + £'~ selection. With this selection, 
87 ZZ — > £ + £~vv candidates are observed in data. 




Figure 3. For t^t'vv candidates in all channels figure (a) shows the axial- E™ 1SB after all selection require- 
ments, except for the axial- E™ 1BB , and figure (b) shows the fractional px difference between E™ 1BB and Px after 
all selection requirements, except for the fractional px difference (the last bin also contains events with frac- 
tional px difference greater than 1). In all plots, the points are data and the stacked histograms show the signal 
prediction from simulation. The shaded band shows the combined statistical and systematic uncertainties. 



4 Signal acceptance 

The Z boson decays to hadrons, neutrinos and charged leptons with branching fractions of 69.9%, 
20.0% and 10.1%, respectively [46]. The two ZZ decay channels considered in this paper, ZZ — > 
£ + £~£' + £'~ and ZZ — > £ + £~~i/v, have branching fractions of 0.45% and 2.69%, respectively 6 , where 
decays involving r leptons are not included in these branching fractions. Some of the ZZ decays 
produce one or more charged leptons which pass through the uninstrumented regions of the detector, 
and as such cannot be reconstructed. In order to measure the total ZZ cross section, the measured 
decays are extrapolated to non-measured parts of the phase-space; this results in the measurement 
being more dependent on theory predictions. Consequently, two types of cross sections are measured: 
fiducial and total. The fiducial cross section is the cross section measured within a restricted phase 
space, and the total cross section is the cross section extrapolated to the total phase space. 

6 The quoted branching fraction to four charged leptons is for the case where both Z bosons are within the mass 
window, so that the 7* contribution can be neglected. 
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Figure 4. (a) Jet multiplicity for the ZZ —> £ + £~ £' + l'~ selection and (b) jet multiplicity for the ZZ — > £ + £~vv 
selection (with all selections applied but the jet veto). The points represent the observed data. In (a) the 
ZZ — > £ + £~ £' + £'~ background is normalized to the data-driven (dd) estimate, while in (b) the histograms show 
the prediction from simulation. The shaded band shows the combined statistical and systematic uncertainty 
on the prediction. 



The total cross section calculation depends on the choice of Z mass range. The cross section is 
calculated using the Z boson natural width rather than the zero-width approximation, and includes 
the mass window requirement (66 to 116 GeV) to remove most of the 7* contamination. The ratio of 
the total cross section calculated with both Z bosons within the mass window to the total cross section 
calculated using the zero- width approximation is 0.953, as the mass window requirement removes some 
of the Z bosons in the tails of the mass distribution. 

4.1 Fiducial region definitions 

The fiducial cross section is restricted to a region which is constructed to closely match the instru- 
mented region and the event selection; for simplicity, only the most inclusive requirements on the 
lepton 77 and are used for the definition of the fiducial phase space. The fiducial cross section cr z A z 
is calculated as: 

zz C zz xC K ' 

which depends on a correction factor given by the number of simulated ZZ^*> events which satisfy the 
full event selection divided by the number of ZZ^ events generated in the fiducial region, Czz\ the 
integrated luminosity, £; the number of selected events, N ^, s ; and the amount of estimated background, 
iVbkg- For the calculation of Czz, final states including pairs of oppositely-charged leptons produced 
from decays of Z — > t + t~ — > t^i~vvvv are included in the number of selected events (numerator) 
since those decays have an identical final state to the signal and are not subtracted as background but 
are excluded from the fiducial region (denominator) because the fiducial regions are defined only with 
ZZ^ decays directly to electrons, muons or neutrinos, depending on the channel. The contribution 
from such r decays is estimated from Monte Carlo simulation to be < 0.1 % for the ZZ — > £ + £~vD 
selection, 0.24±0.01% for the ZZ £+£-£' +£'~ selection and 1.73±0.04% for the ZZ* £ + £-£' + l'- 
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selection. Fiducial requirements are applied at generator level. To reduce the dependence on QED 
radiation, dressed leptons are always used, for which the lepton four-momentum is summed with the 
four- momentum of all photons within AR = 0.1. 

The ZZ — > £ + 1~ ' £' + £'~~ fiducial region is defined using the following requirements: (i) two pairs of 
same-flavour opposite-sign electrons or muons, with each lepton satisfying pi^ > 7 GeV, \rf\ < 3.16 
and at least a distance AR = 0.2 from any other selected lepton, i.e., AR(£i : £2) > 0.2, and (ii) 
both dilcpton invariant masses within the Z mass window. A ZZ* — > £ + £~£' + £'~ fiducial region is 
defined with the same criteria as in the ZZ — > £ + £~£' + £'~ case, except that one dilepton invariant 
mass requirement is relaxed to be greater than 20 GeV. 

The ZZ — > £ + £~ vD fiducial region is defined by requiring: (i) two same- flavour opposite-sign 
electrons or muons, each with > 20 GeV, |?/| < 2.5, with AR(£ + ,£~) > 0.3, (ii) dilepton invariant 
mass close to the Z boson mass: 76 < m e + e - < 106 GeV, (hi) dineutrino invariant mass close to 
the Z boson mass: 66 < m V v < 116 GeV, (iv) no jet with p J T > 25 GeV and \rf\ < 4.5, and (v) 
{\p™ — p^\)/Pt < 0.4 and — p^ ■ p z /p z > 75 GeV. Jets are defined at generator level using the 
same jet algorithm as used in reconstructed events and including all final state particles after parton 
showering and hadronization. 

Fiducial cross sections are calculated using the ZZ £ + £~£ ,+ £'^, ZZ* £+£-£'+£'- and ZZ -» 
£ + £~vv selections, integrated over the corresponding full fiducial phase space volumes. For the ZZ — > 
£ + £~£' + £'~ and ZZ — > £ + £~vv selections the differential fiducial cross sections are derived in bins of 
the leading p 2 ,, A(j>(£ + ,£~) and the mass of the ZZ — > £ + £~£' + £'~ system or the transverse mass of 
the ZZ — !• £ + £~vv system. 

The correction factor, Czz, is determined from Monte Carlo simulations (PowhegBox for the 
ZZ -> £ + £~vv channel and PowhegBox and GG2zz for the ZZ^ -> £ + £~£' + £'~ channel), after 
applying data-driven corrections as described in section 2.1. For the ZZ —> £ + £~£' + £'~ (ZZ* — > 
£+£-£'+£'-) selection it is 0.43 (0.41) for e+e^e+e-, 0.68 (0.69) for fi+ fj,- n + pr and 0.55 (0.53) for 
e + e~/i + /i~ events. For the ZZ — » £ + £~vv selection the correction factor is 0.63 for e + e~vv and 0.76 
for [i + \x~vv events. The correction factors combining all lepton categories within the fiducial region 
are given in table 1 for the three event selections in both decay channels. 



Selection 




Czz 


ZZ -> £+£- 


£'+£'- 


0.552 ± 0.002 ± 0.021 


ZZ* -> £+£ 


-£'+£'- 


0.542 ± 0.002 ± 0.022 


ZZ -+ £+£- 


vv 


0.679 ± 0.004 ± 0.014 



Table 1. Correction factors Czz for each production and decay channel. The first uncertainty is statistical 
while the second is systematic. 



4.2 Extrapolation to the total phase space 

The total ZZ cross section is measured using the ZZ — > £ + £~£' + £'~ and ZZ — > £ + £~vv selections. 
The total cross section is calculated using the fiducial acceptance, Azz (the fraction of ZZ events with 
Z bosons in the Z mass window that fall into the fiducial region) and the branching fraction, BF: 

total _ -V b s — -/Vpkg / . „s 

ZZ ~ A Z z x Czz x £ x BF 1 ' ' 

The fiducial acceptances Azz are estimated from Monte Carlo simulation, using PowhegBox 
for the ZZ — !• £ + £~vv channel and PowhegBox and GG2zz for the ZZ -> £ + £~£' + £'~ channel. The 
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fiducial acceptance of the ZZ — > £ + £ vv channel is much more constrained than the ZZ — > £ + £ £' + £'' 
channel in order to reduce background. Values are given in table 2. 



Selection 






A ZZ 






ZZ -* £+£-£'+£'- 


0.804 


± 


0.001 


± 


0.010 


ZZ -> £+£~vv 


0.081 


± 


0.001 


± 


0.004 



Table 2. Acceptance Azz for the two decay channels used for the measurement of the total ZZ production 
cross section. The first uncertainty is statistical while the second is systematic. 

4.3 Systematic uncertainties 

Table 3 summarizes the systematic uncertainties on Czz and Azz- For Czz in the ZZ^*> — > 
£ + £~£ l+ £'~ selections, the dominant systematic uncertainties arise from the lepton reconstruction 
efficiency, the efficiency of the isolation and impact parameter requirements, and the differences in 
Czz estimated by Sherpa and PowhegBox; uncertainties on the trigger efficiency and the lepton 
energy scale and resolution are small. In the ZZ — > £ + £~vv channel the dominant Czz uncertainties 
are from uncertainties on the lepton reconstruction efficiency, the lepton energy scale and resolution, 
and the missing transverse momentum modelling and jet veto uncertainty; uncertainties on the trigger 
efficiency and due to differences in Czz estimated by Sherpa and PowhegBox also contribute. 

The uncertainties on Czz from the reconstruction efficiency, energy scale and resolution, isolation 
and impact parameter requirements and trigger efficiency are estimated by varying the data-driven 
correction factors applied to simulation by their systematic and statistical uncertainties. The system- 
atic uncertainties on events with extended leptons used in the ZZ^ — > £ + £~£' + £'~ channel are slightly 
higher than in events without them; nevertheless, since their relative contribution is small, the effect 
on the uncertainty of the combined channels is negligible. The generator systematic uncertainty for 
Czz accounts for the effect of choosing a different renormalization and factorization scale and PDF 
set. 

For Azz-, the systematic uncertainties are due to theoretical uncertainties which come from the 
PDFs, the choice of the renormalization and factorization scales, the modelling of the contribution 
from gg initial states and the parton shower model, as given in table 3. For the ZZ — > £ + £~uP channel, 
uncertainties in the efficiency of the jet veto are also taken into account through the calculation of a 
scale factor; the ratio of the jet veto efficiency in data to that in MC simulation is taken from a sample of 
single Z events and then applied to ZZ events [47] . The systematic uncertainties due to the PDFs and 
scales are evaluated with MCFM by taking the difference between the Azz obtained using the CT10 
and MSTW2008 PDF sets, as well as using the 44 CT10 error sets, and by shifting the factorization and 
renormalization scales up and down by a factor of two. An additional uncertainty is assigned to account 
for the effect of different modelling at the generator level. Since the ZZ* — > £ + £~£' + £'~ measurement 
is not used for the total cross section, its Azz acceptance is irrelevant and only uncertainty values 
related to Czz are given. 

The uncertainty on the integrated luminosity is 3.9% [22]. The uncertainty on the background 
estimates is discussed in the following sections. 



- 12 - 



Source 


ZZ -> £+£-£'+£'- 


zz* -> e+t-f+e'- 


ZZ -> l + l~vv 


Czz 


Lcpton efficiency 


3.0% 


3.1% 


1.3% 


Lepton energy /momentum 


0.2% 


0.3% 


1.1% 


Lepton isolation and impact parameter 


1.9% 


2.0% 


0.6% 


Jct+£^ liss modelling 






0.8% 


Jet veto 






0.9% 


Trigger efficiency 


0.2% 


0.2% 


0.4% 


PDF and scale 


1.6% 


1.5% 


0.4% 


Azz 


Jet veto 






2.3% 


PDF and scale 


0.6% 




1.9% 


Generator modelling and parton shower 


1.1% 




4.6% 



Table 3. Summary of systematic uncertainties, as relative percentages of the correction factor Czz or the 
acceptance of the fiducial region Azz- Dashes indicate uncertainties which are not relevant. 



5 Background estimation 

5.1 ZZ« l+trl'+e- background 

Background to the ZZ^ — > £ + £~f' + £'~ signal originates from events with a Z (or W) boson decaying 
to leptons accompanied by additional jets or photons (W/Z + X), from top-quark production and 
from other diboson final states. Such events may contain electrons or muons from the decay of heavy- 
flavoured hadrons, muons from in-flight decay of pions and kaons, or jets and photons misidentified as 
electrons. The majority of these background leptons are rejected by the isolation requirements. 

The background estimate follows a data-driven method in which a sample of events containing 
three leptons satisfying all selection criteria plus one 'lepton-likc jet' is identified; such events are 
denoted as £££j. For muons, the lepton-likc jets are muon candidates that fail the; isolation requirement 
or fail the impact parameter requirement but not both. For electrons with |?7| < 2.47, the lepton-like 
jets are clusters in the electromagnetic calorimeter matched to inner detector tracks that fail either 
the full electron selection or the isolation requirement but not both. For electrons with |?7| > 2.5, 
the lepton-like jets are electromagnetic clusters that are reconstructed as electrons but fail the tight 
identification requirements. The events are otherwise required to satisfy the full event selection, 
treating the lepton-like jet as if it were a fully identified lepton. The background is then estimated by 
weighting the £££j events by a measured factor /, which is the ratio of the probability for a non-lepton 
to satisfy the full lepton selection criteria to the probability of satisfying the lepton-like jet criteria. 
The background in which two selected leptons originate from jets is treated similarly, by identifying 
a data sample with two leptons and two lepton-like jets; such events are denoted as ££jj. The total 
number of expected background £ + £~£' + £'~ events, iV(BG), is calculated as: 

7V(BG) = [N(£££j) - N(ZZ)} x / - N(Ujj) x f (5.1) 

where double counting from £££j and ££jj events is accounted for, and the term N(ZZ) is a Monte 
Carlo estimate correcting for contributions from signal ZZ^ — > £ + £~£ /+ £'~ events having a real lepton 
that is classified as a lepton-like jet (the equivalent correction to the term N(££jj) is negligible). 
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i i 
e e e T e 


1 1 

A 4 A 4 At A* 


i i 
e e fi n 


i + i~m + d - 


(+) N(£££j) x / 


1.63 ± 0.34 


0.21 ± 0.21 


1.84 ± 0.40 


3.67 ±0.57 


(-)N(ZZ) x / 


0.17 ±0.13 


12+ - 20 


0.34 ±0.21 


0.63 ±0.32 


(-)N(Ujj) x f 


0.96 ±0.10 


0.33 ±0.16 


0.83 ±0.09 


2.12 ±0.21 


Background estimate, 7V(BG) 


0.5lg;|(stat.) 


< 0.64 


0.7 ± 0.7(stat.) 


0.9lJ;J(stat.) 




±0.3(syst.) 




±0.6 (syst.) 


±0.7(syst.) 



Table 4. Expected number of background events for the ZZ — > £ + £~ £' + 1'~ selection in 4.6 fb _1 of data, for 
the individual decay modes (columns 2, 3 and 4) and for their combination (last column). If the central value 
of the estimate is negative, the upper bound on the number of events in that channel is derived as detailed in 
section 5.1. 

The factor / is measured in a sample of data selected with single-lepton triggers which contain a 
Z boson candidate: a pair of isolated same-flavour opposite-sign electrons or muons. In these selected 
events, / is measured, using the lepton and lepton-like jet candidates not assigned to the Z boson, 
as the ratio of the number of selected leptons to the number of lepton-like jets, after correcting for 
expected true lepton contributions from WZ and ZZ events using simulation. Independent values as 
a function of the 77 and of the lepton-like jet are measured, which are then combined assuming they 
are uncorrelated. The factor / is found to vary from 0.33 ± 0.01 (0.26 ± 0.02) below pt = 10 GeV 
to 0.09 ± 0.02 (0.46 ± 0.20) above pt = 50 GeV for electrons (muons). The quoted uncertainties are 
statistical. Then, with the same procedure, a value for / is also derived using the simulated samples 
of background processes. The difference between the value of / derived in data and in simulation is 
assigned as a systematic uncertainty on /. The statistical and systematic uncertainties are then added 
in quadrature to derive a combined uncertainty on /, which varies as a function of from 14% (19%) 
below 10 GeV to 22% (51%) above 50 GeV for electrons (muons). For the muons, the total uncertainty 
on / is dominated by its statistical uncertainty. The background estimates for the ZZ — > £ + £~£' + £'~ 
and ZZ* ->■ £+£-£'+£'- selections are 0.9±J:| (stat.) ±0.7(syst.) and 9.1 ± 2.3(stat.) ± 1.3(syst.) events, 
respectively, as shown in tables 4 and 5. The statistical uncertainty on the background estimate comes 
from the statistical uncertainty on the numbers of £££j, ££jj and ZZ^*' — > £ + £~ events used in 
eq. 5.1. The systematic uncertainty results from the combined uncertainty on /. In cases where the 
overall estimate is negative, the background estimate is described using a truncated Gaussian with 
mean at zero and standard deviation equal to the estimated statistical and systematic uncertainties 
added in quadrature. 

The extra background induced by the use of the extended leptons in the ZZ^ — > £ + £~£' + £'~ 
channel is estimated to be negligible in the ZZ — > £ + £~~ £' + £'~ selection, and about 20% (2 events out 
of the 9.1 estimated, compared to a signal gain of about 10.6 events out of the 64.4 expected) in the 
ZZ* £+£-£'+£'- selection. 

The background is also estimated purely from the simulated samples of background processes, and 
is predicted to be 1.5 ± 0.4 events for the ZZ — > £ + £~£' + £'~ selection and 8.3 ± 1.3 events for the 
ZZ* — > £ + £~£' + £'~ selection, with uncertainties being statistical only. These estimates compare well 
with the data-driven results given in tables 4 and 5. According to the estimate from simulation, the 
dominant source of background is Z+jets events, with only about a 10% to 20% contribution from 
other diboson channels (WZ and WW), and a negligible contribution from events with top quarks. 

Differential background distributions are determined by first deriving the shape of the distributions 
from the background MC samples. This is achieved by selecting events where one Z candidate is 
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i i 
e e e T e 


1 1 

A 4 A 4 At A* 


1 1 


£+£-£'+£' ~ 


(+) N(£££j) x / 


8.85 ± 0.98 


0.21 ± 0.21 


10.63 ± 1.06 


19.70 ± 1.46 


(-)N(ZZ) x / 


0.29 ±0.18 


n on+ - 25 

u - zu -0.20 


0.56 ±0.28 


1.05 ±0.42 


(-)N(Ujj) x p 


4.24 ±0.23 


1.10 ±0.31 


4.24 ±0.23 


9.58 ±0.45 


Background estimate, iV(BG) 


4.3 ± 1.4(stat.) 


< 0.91 


5.8 ± 1.6(stat.) 


9.1 ± 2.3(stat.) 




±0.6(syst.) 




±0.9 (syst.) 


±1.3(syst.) 



Table 5. Expected number of background events for the ZZ* — > £ + £~ £' + £'~ selection in 4.6 fb _1 of data, for 
the individual decay modes (columns 2, 3 and 4) and for their combination (last column). If the central value 
of the estimate is negative, the upper bound on the number of events in that channel is derived as detailed in 
section 5.1. 



Process 




e + e~E™ iss 




£+£-E™ iss 


ti, Wt, WW, Z - 


-> T+T 


8.5 ±2.1 ±0.5 


10.6 ±2.6 ±0.6 


19.1 ± 2. 3± 1.0 


WZ 




8.9 ±0.5 ±0.4 


11.9 ±0.5 ±0.3 


20.8±0.7±0.5 


Z — » e + e~ 


±jets 


2.6 ±0.7 ±1.0 


2.7±0.8± 1.2 


5.3 ±1.1 ±1.6 


W+ jets 




0.7 ±0.3 ±0.3 


0.7 ±0.2 ±0.2 


1.5 ± 0.4 ± 0.4 


W~/ 




0.1 ±0.1 ±0.0 


0.2 ±0.1 ±0.0 


0.3 ±0.1 ±0.0 


Total 




20.8 ±2.3 ±1.2 


26.1 ± 2.8 ± 1.4 


46.9 ± 4.8 ± 1.9 



Table 6. Expected number of background events to the ZZ — > £ + £~vi> channel in 4.6 fb _1 of data, for the 
individual decay modes (columns 2 and 3) and for their combination (last column). The first uncertainty is 
statistical while the second is systematic. 

required to satisfy the nominal lepton selection, while the other Z candidate is formed by leptons 
satisfying relaxed criteria for the isolation requirements and transverse impact parameter significance. 
The shape determined in this way is then scaled such that the total number of events in the distribution 
is equal to the data-driven background estimate shown in tables 4 and 5. 

5.2 ZZ £+£-vv background 

There are several sources of background to the ZZ —> £ + £rvv channel. Processes such as ti, WW, Wt 
or Z — > t + t~ production give two true isolated leptons with missing transverse momentum. Diboson 
W Z events in which both bosons decay leptonically have three charged leptons, but if one lepton from 
a W or Z boson decay is not identified, the event has the same signature as the signal. Production of 
a Z boson in association with jets gives two isolated leptons from the Z boson decay and may have 
missing transverse momentum if the jet momenta are mismeasured. Finally, production of a W boson 
in association with jets or photons may satisfy the selection requirements when one of the jets or 
photons is misidentified as an isolated lepton. All of the backgrounds are measured with data-driven 
techniques except for WZ and Wj. The total background is estimated to be 46.9 ± 4.8 ± 1.9 events 
as summarized in table 6. 

5.2.1 Backgrounds from it, Wt, WW and Z t+t~ 

The contributions from ti, Wt, WW and Z — » t + t~ processes are measured by extrapolating from 
a control sample of events with one electron and one muon (instead of two electrons or two muons) , 
which otherwise satisfy the full ZZ — > £ + £~vv selection. This sample is free from signal events. 
The extrapolation from the e/z channel to the ee or \x\x channel uses the relative branching fractions 
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(2 : 1 : 1 for e/i : ee : /x/x) as well as the ratio of the efficiencies e ee or of the ee or /x/x selections 
to the efficiency e efl of the efi selection, which differs from unity due to differences in the electron and 
muon efficiencies. 

For the electron channel, this is represented by the equation: 

^V e b e ks = « ata - KT) x \ x e — (5.2) 

» eM is the number of observed e[i events and JV|™ is the number expected events from 
processes other than tt, Wt, WW and Z -> t+t~ (WZ, ZZ, VK+jct, Z+jet and W/~f). Therefore, 
(Ndata _ 7V|™) is the estimate of tt, Wt, WW and Z -> t+t~ production in the control sample. The 
efficiency correction factor, e ee /e e ^, corrects for the difference between electron and muon efficiency. 
The efficiency correction factor is measured in data using reconstructed Z — > P + l~ events, as 



/ydata 

(5.3) 



where iV e d e ata and N^ ta are the number of observed ee or /x/x events in the Z boson mass window, 
respectively, after all lepton selection requirements and the Z boson mass window requirement are 
applied. A parallel argument gives Nj^ s . This procedure is repeated in bins of in order to obtain 
the pt distribution of the tt, Wt, WW and Z — > t+t~ background. 

The dominant uncertainty is statistical (25%), due to the limited number of events in the control 
samples. Additional uncertainties are due to systematic uncertainties in the normalization of the 
simulated samples used to correct the e/x contribution (5.5%) and the systematic uncertainty in the 
efficiency correction factor (4.5%). 

5.2.2 Background from WZ production with leptonic decays 

Events from leptonic WZ decays may result in an £ + £~ E™ ss signature when one lepton from the W 
or Z boson is not reconstructed. The contribution from this process is estimated using the simulated 
samples described in section 2.1. The estimate is checked using a control region with three high- 
Pt isolated leptons. The two dominant processes that contribute to this control region are WZ 
and Z+jets production, where the WZ boson pair decays to three leptons and a neutrino and the 
-Z+jets contribution has two real leptons from the Z decay and a misidentified lepton from the jet. 
The technique used to estimate the background in the ZZ^ — > i + (r £' + l'~ channel is also used to 
normalize the contribution from Z+jets in the three-lepton control region. The WZ Monte Carlo 
expectation is consistent with the data. The systematic uncertainties are estimated in the same way 
as for signal Monte Carlo events. 

5.2.3 Background from Z bosons with associated jets 

Occasionally events with a Z boson produced in association with jets may have large amounts of 
missing transverse momentum due to mismeasurement of the momenta of the jets. This background 
is estimated using events with a high-p T photon and jets as a template, since the mechanism for 
large missing transverse momentum is the same as in Z+jets events. The events are reweighted such 
that the photon Et matches the observed Z boson pt and are normalized to the observed Z + jets 
yield. The procedure is repeated in bins of p^ in order to obtain the pt distribution of the Z+jets 
backgrounds. The largest systematic uncertainty is due to the subtraction of W-y, Zj, tt and W — > ev 
contributions to the 7+jets sample, which is 33% in the ee channel and 37% in the /x/x channel. 
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5.2.4 Background from events with a misidentified lepton 

A small contribution to the selected sample is due to events in which one of the two leptons comes 
from the decay of a W or Z boson (called 'real' below) and the second is a 'fake', corresponding both 
non-prompt leptons and misidentified ir° mesons or conversions. 

The dominant fake-muon mechanism is the decay of heavy-flavoured hadrons, in which a muon 
survives the isolation requirements. In the case of electrons, the three mechanisms are heavy-flavour 
hadron decay, light-flavour jets with a leading 7r° overlapping with a charged particle, and conversion 
of photons. Processes that contribute are top-quark pair production, production of W bosons in 
association with jets and multi-jet production. 

The 'matrix method' [48] is applied to estimate the fraction of events in the signal regions that 
contain at least one fake lepton. The method measures the number of fake leptons in background- 
dominated control regions and extrapolates to the ZZ selection region using factors measured in data. 
The shape of the background is provided by taking the background as uniformly distributed among 
the bins and treating each bin as statistically uncorrelated. The dominant systematic uncertainty is 
due to the uncertainty on the extrapolation factors and the limited numbers of events in the control 
samples, giving a total uncertainty of 63% and 44% in the ee and channels, respectively. 

6 Results 

Three types of measurements are presented: 

• integrated fiducial and total ZZ cross sections; 

• differential cross sections normalized to the overall measured cross sections for the and 
A(f)(£ + of the leading Z boson, and the mass (transverse mass 7 ) of the ZZ system for 
the ZZ -> £+£-£' +f" (ZZ ->• £+£~vv) selection; and 

• limits on the anomalous nTGCs. 
6.1 Cross section measurements 

The expected and observed event yields after applying all selection criteria are shown in table 7 for 
both channels. Figure 4 shows the jet multiplicity in selected ZZ — > £ + £~£' + £'~ and ZZ — > £ + £~vv 
events before the jet veto is applied. Figures 5 and 6 show the transverse momentum and mass of 
the ZZ system in selected ZZ £+£-£'+£'- and ZZ* £+£-£'+£'" events respectively. Figure 7 
shows the transverse momentum and mass of the two-charged-lepton system in selected ZZ — > £ + £~ vv 
events. The shapes of the distributions are consistent with the predictions from the simulation. 

The ZZ(*' — > £ + £~£' + £'~ and ZZ —> £ + £~vv fiducial cross sections are determined using a maxi- 
mum likelihood fitting method, taking into account the integrated luminosity and the Czz correction 
factors discussed in section 4. A Poisson probability function is used to model the number of ex- 
pected events, multiplied by Gaussian distribution functions which model the nuisance parameters 
representing systematic uncertainties. The measured fiducial cross sections are: 




(syst.) ± 1.0 (hum.) fb, 



(syst.) ± 1.2 (lumi.) fb, 
(syst.) ± 0.5 (lumi.) fb. 




^{m z ) 2 + (pf ) 2 + y'K) 2 + (E^)A - (pf + Hj? iss ) 



2 
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ZZ K > -> r£ t + i 


4- — _L — 

e e e e 


4- — -4- — 

MM MM 


-i_ _ -i_ _ 
e T e H (J, 


fi-\- n — nf A- fit — 

r£ t ^l 


Observed ZZ 

s~\ i l r~7 rv-ii 

Observed ZZ 


16 
21 


23 
30 


27 
33 


66 
84 


Expected ZZ signal 
Expected ZZ* signal 


10.3 ± 0.1 ± 1.0 
12.3 ± 0.2 ± 1.2 


16.5 ± 0.2 ± 0.9 
20.5 ± 0.2 ± 1.1 


26.7 ± 0.2 ± 1.7 
31.6 ± 0.3 ± 2.0 


53.4 ± 0.3 ± 3.2 
64.4 ± 0.4 ± 4.0 


Expected ZZ background 
Expected ZZ* background 


0.5 ± 0.6 ± 0.3 
4.3 ± 1.4 ± 0.6 


< 0.6 

< 0.9 


0.7 ± 0.7 ± 0.6 
5.8 ± 1.6 ± 0.9 


0.9 ± 1.1 ± 0.7 
9.1 ± 2.3 ± 1.3 




ZZ -> g+l-vv 


e+e-E™ iss 








Observed ZZ 


35 


52 




87 


Expected ZZ signal 


17.8 ± 0.3 ± 1.7 


21.6 ± 0.3 ± 2.0 




39.3 ±0.4 ±3.7 


Expected ZZ background 


20.8 ±2.3 ±1.2 


26.1 ± 2.8 ± 1.4 




46.9 ± 4.8 ± 1.9 



Table 7. Summary of observed ZZ -> £ + T £'+£'-, ZZ* -> £+£~ £'+£'- and ZZ -J- £ + £~w candidates in the 
data, total background estimates and expected signal for the individual decay modes (columns 2 to 4) and 
for their combination (last column). The quoted uncertainties and limits represent 68% confidence intervals; 
the first uncertainty is statistical while the second is systematic. The uncertainty on the integrated luminosity 
(3.9%) is not included. 
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Figure 5. (a) Transverse momentum p^ z and (b) invariant mass 



of the four-lepton system for the ZZ 



selection. The points represent the observed data and the histograms show the prediction from simulation, 
where the background is normalized to the data-driven (dd) estimate as described in section 5.1. The shaded 
band shows the combined statistical and systematic uncertainty on the prediction. 



where £ + £ £' + £' refers to the sum of the e + e e + e , e + e li + \i and final states and 

£ + £~vv refers to the sum of the e + e~E™ ms and (i + n~E!f lss final states 8 . The expected SM fiducial 

8 The ZZ — > £+l~vi> fiducial region is more restricted compared to the ZZ^ — > £~ £'^~ £'~ channel. 
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Figure 6. (a) Transverse momentum pf z and (b) invariant mass m zz of the four-lepton system for the ZZ* 
selection. The points represent the observed data and the histograms show the prediction from simulation, 
where the background is normalized to the data-driven (dd) estimate. The shaded band shows the combined 
statistical and systematic uncertainty on the prediction. 
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Figure 7. (a) Transverse momentum p z and (b) mass mz of the two-charged-lepton system for the ZZ — >• 
l + l~vv selection. The points represent the observed data and the histograms show the prediction from 
simulation. The shaded band shows the combined statistical and systematic uncertainty on the prediction. 



cross sections, derived from PowhegBox and GG2zz, are: 



ofz^i+ri'+l'- = 20 - 9 ± °- 1 ( stat -) -0:9 (theory) fb, 
°f?™ t+f - t , +f ,- = 25.6 ± 0.1 (stat.) t\i (theory) fb, 



J "' SM - 12.5 ± 0.1 (stat.) t\i (theory) fb. 



The measured cross sections are compatible with these theoretical values. 
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The total ZZ cross section is calculated by extrapolating to the full phase space while each Z boson 
is required to have a mass within the Z mass window. Both ZZ — > £ + £~£' + £'~ and ZZ — > £ + £~vv 
events are combined in the maximum likelihood fit, taking into account the known Z branching 
fractions [46] and the Azz kinematic and geometrical acceptances (section 4). Correlated systematic 
uncertainties between the ZZ — > £ + £~£' + £'~ and ZZ — > £ + £~vv channels are taken into account in 
the fit using a single Gaussian for the nuisance parameter for each source of correlated uncertainty. 
The measured value of the total ZZ cross section is: 

afz = 6.7 ± 0.7 (stat.) ±g | (syst.) ± 0.3 (lumi.) pb. 

The result is consistent within errors with the NLO Standard Model total cross section for this process 
of 5.89ito'i8 pb, where the quoted theoretical uncertainties result from varying the factorization and 
rcnormalization scales simultaneously by a factor of two and from using the full CT10 PDF error set. 

6.2 Differential cross sections 

The differential cross sections present a more detailed comparison of theory to measurement, allowing 
a generic comparison of the kinematic distributions to new theories. Variables which are sensitive to 
new phenomena, such as p z , m zz and A<fi(£ + ,£~), are used with bin boundaries chosen to maximize 
sensitivity to nTGCs. At the same time, the bin widths were chosen to be commensurate with the 
resolution. 

The measured distributions are unfolded back to the underlying distributions, accounting for the 
effect of detector resolution, efficiency and acceptance, within the fiducial region of each measurement. 
The unfolding procedure is based on a Bayesian iterative algorithm [49]. The algorithm takes as input 
a prior for the kinematic distribution and iterates using the posterior distribution as prior for the 
next iteration. The initial prior is taken from the signal Monte Carlo expectation calculated using 
the PowhegBox generator and three iterations are performed. The uncertainty on the unfolded 
distributions is dominated by the statistical uncertainty, which is about 30% in most bins. The 
systematic uncertainty is no more than 5% in any bin. The dependence of the unfolded cross sections 
on the choice of the initial prior is tested by unfolding the measured distributions using a different 
generator (Sherpa). The difference between the two is taken as a systematic uncertainty to account 
for differences in generator modelling (e.g. QCD radiation). The difference in unfolded distributions 
between three iterations and four iterations is much lower than the statistical uncertainty and it is 
taken as a further uncertainty on the unfolding procedure. Systematic uncertainties related to detector 
effects (e.g. lepton reconstruction efficiency) are evaluated using pseudo-experiments. 

Figures 8 to 10 show the differential cross sections normalized to the fiducial cross sections for the 
p z and A(f>(£ + ,£~) of the leading Z boson, and for the mass (transverse mass) of the ZZ system for 
the ZZ £+£-£'+£'- (ZZ £+£^vD) selection. The Standard Model prediction is consistent with 
the measurement in each case. 

6.3 Anomalous neutral triple gauge couplings 

Anomalous nTGCs for on-shcll ZZ production can be parameterized by two CP-violating (fX) and two 
CP-conserving (f§) complex parameters (where V = Z, 7) which are zero in the Standard Model [3]. 
A form-factor parameterization is introduced leading to couplings which vanish at high parton centre- 
of-mass energy y/§: fV — /^/(l + s/A 2 ) n , ensuring partial-wave unitarity. Here, A is the energy scale 
at which physics beyond the Standard Model would be directly observable, are the low-energy 
approximations of the couplings, and n is the form-factor power. Values of n — 3 and A = 3 TeV are 
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Figure 8. Unfolded ZZ fiducial cross sections in bins of the pr of the leading Z boson for (a) the ZZ — > 
£ + £~£' + £'~ selection, where a discontinuity is indicated by the parallel pairs of lines, and (b) the ZZ — > £ + £~vv 
selection. 
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Figure 9. Unfolded ZZ fiducial cross sections in bins of the A(j>(£ + ,£ ) of the leading Z boson for (a) the 
ZZ -> £+£-£'+£'- selection and (b) the ZZ -> £ + £~vu selection. 



chosen, so that expected limits are within the values allowed by requiring that unitarity is not violated 
at LHC energies [3]. The results with an energy cutoff A = oo (i.e. without a form factor) are also 
presented as a comparison in the unitarity violating scheme. 

Limits on anomalous nTGCs are determined using the observed and expected numbers of ZZ — > 
£ + £~£' + £'~ and ZZ — > l + l~vv events binned 9 in p^, as seen in table 8. Figure 11 shows the observed 
p|i distributions, together with the SM expectation and the predicted distributions for nTGC values 
close to the previous limits obtained by ATLAS [13]. Using an increased data sample compared 

9 The raw (i.e. not unfolded) differential event yields are used, to avoid introducing theory dependence. 
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Figure 10. Unfolded ZZ fiducial cross sections in bins of (a) m for the ZZ -> t + ri' + l'~ selection and 
(b) mf z for the ZZ -» l + t~vv selection. 





Expected background 


Expected ZZ signal 


Observed events 


zz -> e+r e +f- 


< pf < 60 GeV 


0.6 ±0.8 ±0.5 


27.9 ± 0.2 ± 2.0 


28 


60 < pf < 100 GeV 


0.2 ±0.2 ±0.2 


14.6 ± 0.2 ± 1.2 


25 


100 < pf < 200 GeV 


0.1 ±0.1 ±0.1 


9.3 ± 0.1 ± 0.9 


11 


pf > 200 GeV 


0.01 ±0.01 ±0.01 


1.6 ± 0.1 ± 0.3 


2 


zz -> t+trvv 


50 < p| < 90 GeV 


26.0 ± 4.5 ± 1.1 


13.6 ±0.2 ±1.3 


42 


90 < pf < 130 GeV 


16.0 ±2.8 ±0.7 


15.7 ±0.3 ±1.7 


29 


pf > 130 GeV 


4.9 ± 1.8 ±0.2 


10.1 ± 0.1 ± 1.5 


16 



Table 8. Total background, expected signal and observed events as a function of the pt of the leading Z for 
the ZZ —> l + l~ l' + £'~ and ZZ —> l + l~uv selections. For the expected signal and background events, the first 
uncertainty is statistical and the second is systematic. 



with our previous measurement, including the ZZ — > £ + £^vD channel, and exploiting the differential 
event yields, the precision is expected to improve by about a factor of five. The dependency of the 
couplings on the expected number of events in each p^ bin is parameterized using fully simulated 
events, generated with Sherpa [26], subsequently reweighted using the Baur-Rainwater [3, 50] and 
BHO [51] MC generators. The next-to-leading-order matrix elements with their nTGC dependence 
have been extracted from the BHO MC generator for 2^5 events and the Baur-Rainwater MC 
generator for 2 — !• 4 events and introduced into a framework [52] that enables a calculation of the 
amplitude given the four vectors and the identity of the incoming and outgoing particles from the 
hard process. 

Confidence intervals for the anomalous triple gauge couplings are determined using the maximum 
profile likelihood ratio. Limits are set on each coupling, assuming all of the other couplings are zero 
(as in the Standard Model), and on pairs of couplings assuming the remaining two couplings are zero. 
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Figure 11. The leading Z boson transverse momentum distributions for (a) the ZZ — > £ + 1~ £' + £'~ selection 
and (b) the ZZ — > £ + £~uv selection. The observed distributions are shown as filled circles, the SM expected 
signal and background are shown as filled histograms, and the predicted distributions for four different nTGC 
samples with form factor scales of A = 3 TeV and nTGC coupling values set near the edge of the exclusion set 
in the 1 fb _1 analysis [13] are shown as dashed lines. 



A 


J 40 


J 40 


ti 

J 50 


J 50 


3 TeV 


[-0.022,0.023] 


[-0.019,0.019] 


[-0.023,0.023] 


[-0.020,0.019] 


oo 


[-0.015,0.015] 


[-0.013,0.013] 


[-0.016,0.015] 


[-0.013,0.013] 



Table 9. One-dimensional 95% confidence intervals for anomalous neutral gauge boson couplings, where the 
limit for each coupling assumes the other couplings are fixed at zero, their SM value. Limits are presented 
for form factor scales of A = 3 TeV and A = oo and include both statistical and systematic uncertainties; the 
statistical uncertainties are dominant. 



The profile likelihood ratio is calculated for the data, and also for 10000 pseudo-experiments generated 
using the expected number of events at each point in the one- or two-dimensional nTGC parameter 
space. A point is rejected if more than 95% of the pseudo-experiments have a larger profile likelihood 
ratio value than the one observed in data. The systematic errors are included as nuisance parameters. 

The resulting limits for each coupling are listed in table 9. Two-dimensional 95% confidence 
intervals 10 are shown in figure 12. The one-dimensional limits are more stringent than those derived 
from measurements at LEP [8] and the Tevatron [9] and previously by ATLAS [13]; it should be noted 
that the limits from LEP do not use a form factor, and those from the Tevatron use A = 1.2 TeV. A 
comparison of the LHC limits with those derived from LEP and Tevatron is shown in figure 13. 

7 Conclusions 

A measurement of the ZZ^ production cross section in LHC proton-proton collisions at ^/s = 7 
TeV is presented with data collected by the ATLAS detector, using the ZZ ( * } -> £+£-£'+£'- and 

10 Since most of the sensitivity of the measurement is contained in a single bin, the likelihood ratio used to obtain the 
two-dimensional limits has one effective degree of freedom. 
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ZZ — > £ + l~vv decay channels. Fiducial cross sections are measured for three production and decay 
selections, and the results are compatible with the SM expected cross sections. Using the ZZ — > 
£ + t~£' + £'~ and ZZ — > £ + £~vv selections, the total ZZ production cross section is determined to be: 

crU = 6.7 ± 0.7 (stat.) ±g:| (syst.) ± 0.3 (lumi.) pb. 

The result is statistically consistent with the NLO Standard Model prediction of 5.891° 'ig pb, calcu- 
lated with Z bosons with a mass between 66 and 116 GeV, and supersedes the previous measurements 
made with part of the same dataset [13]. Unfolded distributions of the fiducial cross sections are 
derived for the and A(j)(£ + ,£-) of the leading Z boson and for m zz in the ZZ £ + £-£'+£'- 
selection and the in the ZZ — > l + £~vi> selection. 

The event yields as a function of the px of the leading Z boson for the ZZ — > £ + l~£' + £'~ and 
ZZ — > £ + £~vD selections are used to derive 95% confidence intervals for anomalous neutral triple gauge 
boson couplings. These limits are more stringent than those derived from measurements at LEP [8] 
and the Tevatron [9]. They improve the previous published results from ATLAS [13] by approximately 
a factor of five and supersede them. 
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Figure 12. Two-dimensional triple gauge coupling limits for form factor scale A = oo. The one-dimensional 
triple gauge coupling limits are shown as vertical and horizontal lines inside the two-dimensional ellipses, whose 
shape is determined by the theoretical correlations. For each two-dimensional limit the other TGC parameters 
are assumed to be zero. Since most of the sensitivity of the measurement is contained in a single bin, the 
likelihood ratio used to obtain the two-dimensional limits has one effective degree of freedom. 
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Figure 13. Anomalous nTGC 95% confidence intervals from ATLAS, LEP [8] and Tevatron [9] experiments. 
Luminosities, centre-of-mass energies and cut-offs A for each experiment are shown. 
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